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1 Introduction 



Hydrodynamics is the natural framework for the nuclear matter at finite temperature and 
at long distance scales. A novel nuclear matter state has been discovered in the Relativistic 
Heavy Ion Collision (RHIC) experiment [1] (for review see for example [2]). It is widely 
believed that a deconfinement phase transition happens and the corresponding deconfined 
state is the quark gluon plasma (QGP) [3] (for early work prior to QCD, see [4]). It has been 
estimated that the critical density and critical temperature are nc ~ IGeV/ fm 3 , ten times 
the density of the normal nuclear matter, and Tq ~ 170MeV at which a s ~ 1. So the relia- 
bility of the perturbative description of QGP is highly suspicious. To apply hydrodynamics 
in this case requires the consideration of the propagation of the color charges. Therefore 
the corresponding liquid is "colored". More interestingly, experiment shows that the shear 
viscosity of the deconfined nuclear matter is close to the lower bound set by theory based on 
AdS/CFT [5]. This implies that a perfect fluid is a good approximation to the real QGP. 
The framework for a perfect color liquid has been developed in [6] (for review see [7]). 

In using the nonabelian fluid mechanics to explore the physics of QGP, the natural first 
step is to discuss the classical solution space. So we ask whether there exists any topological 
nontrivial solution. In fact within the framework of the nonabelian fluid, the configuration 
of fluid is described by a field g which takes values in the color group, say SU(n). Because 
QGP only exists in finite space region, a natural boundary condition is that g should go to a 
constant at spacial infinity. Therefore, the fluid configurations are classified by Ii^,{SU{n)) = 
Z. Moreover, in each topological class, we hope to check whether there exists a configuration 
minimizing the total energy. Such solutions are referred to as skyrmions related to color 
group or simply color skyrmions. In our previous paper [8], the existence of color skyrmion is 
shown for a particular choice of the hydrodynamic Hamiltonian. In this paper, we continue the 
discussion of the existence of color skyrmions in both time-independent and time-dependent 
cases. 

We will consider the nonabelian fluid system whose equation of state (EOS) follows the 
so-called j-law, namely the pressure density is proportional to the energy density: 



From the relativistic fluid mechanics, we know that the case 7 = 4/3 corresponds to radiation 
(p = e/3) while the case 7 = 1 corresponds to dust (p = 0). The up-to-date results on EOS 
for the isentropic expansion of QGP come from lattice simulation [9]. In the temperature 
region relevant to RHIC, for example, a phenomenological parametrization of the resulting 
equation of state is given by 



P = (7 - 



(1) 




(2) 
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We see it is a smooth transition from radiation to dust. In addition, we will take 7 = 2 as a 
heuristic example. Consequently, we will restrict 1 < 7 < 2 in this paper. Our major results 
can be summarized as follows. For the time-independent case, 

• 7 = 2: we showed the existence of color skyrmion for this case in a particular ansatz for 
configuration in previous work. In this paper, the existence is shown for more general 
cases. 

• 6/5 < 7 < 5/3: the topological configuration are in general unstable. The total energy 
gets minimized by infinite dilution in the fluid and the final value of energy is zero. 

• 7 = 1: the energy is minimized to a finite value when the fluid is infinitely diluted. 

• l<7<6/5 and 5/3 < 7 < 2: unstable topological configurations exist; however, the 
existence of any metastable soliton cannot be ruled out. 

For time-dependent case, both oscillating and expanding solutions are found. 

The organization of this paper is the following. In Sect. 2, preliminaries materials on 
nonabelian fluid mechanics are reviewed. In Sect. 3 and 4, static color skyrmions are thor- 
oughly considered for different 7. In Sect. 5, time-dependent solutions are given. In Sect. 6, 
conclusions are given and some open issues are addressed. 



2 Preliminaries 

We describe the internal degrees of freedom for a single colored particle by a group element 
g in the color group SU (n) . And the corresponding phase space structure is determined by 
the symplectic potential: 

6 = -ipsTriTsg^dg) (3) 
from which a key Poisson bracket is determined 

{ Ps ,g} = -igT s . (4) 

T s is a generator in Cartan subalgebra and p s is the conjugate momentum. The index s runs 
over all the generators T s . The dynamics of color degrees of freedom is therefore determined 
by the Lagrangian 

L = -ip s Tr(T s g~ l g) - H(p s ,g) (5) 

where H is the Hamiltonian. 
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The generalization to many-particle system gives the nonabelian fluid Lagrangian density: 

£ = X>?ov-F({n 8 }) (6) 

s 

where j£ is the color current, = jsjsn and F an invariant potential function. And 

«V = -iTr(T a g-%g) (7) 

is the generalized velocity field. The action in a general metric background is given by 
S = f d^Xy/gC. In this paper, we will only consider rank-one group SU(2); therefore, the 
index s will be omitted hereafter. 

Recall that the energy-momentum tensor of an ideal fluid is of the standard form 

r" = ( £ + p)//-p/" (8) 

where is the unit velocity field such that u^u^ = 1. Now we identify in the following 
way: 

f = nu^. (9) 

Note an on-shell condition is 

nco^ = F'j^. (10) 
With the identity Sg = gg^ u 5g^ v and the definition of the energy-momentum tensor, 

2 5S 

T ^ = ii 

V9og^ 

it can be verified that T^ v for the nonabelian fluid system in Eq. (6) is in accordance with 
the form for an ideal fluid in Eq. (8) and that 

e = F, p = nF' — F. (12) 

So we claim that the Lagrangian in Eq. (6) does describe an ideal fluid system. Eq. (12) 
shows that for the 7-law of Eq. (1), we solve 

F{n) = -n\ (13) 

7 

with a a dimensional constant. 

Remember the physical boundary condition for g as a field defined in space M. 3 in the 
previous introductory section. In accordance with the boundary condition, the configuration 
space is classified by the mapping class g : S 3 — > SU(n). We do not need to specify the 
topological number in this paper. The assumption that the configuration carries certain 
topological charges is enough in practice below. 
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3 Variational Analysis of Time-Independent Configurations 



We discuss the existence of the static soliton solution with defined topological charge in the 
nonabelian fluid in this section. 



3.1 General Setting 

To do so, it is convenient to use Hamiltonian. In fact, the Lagrangian in Eq. (6) can be 
reexpressed as C = pujQ — 7i where p = j°. The Hamiltonian density is given by 



H=]-u + F (14) 



and there is an on-shell constraint: 



6H SF F' -> 

= UJ+ =UJ j = 0. (15) 

5j Sj n 



With some simple algebra, we get a dynamical system 
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n=^ + F, p 2 = n 2 + !^ 2 (16) 

where to = \co\. In principle, the invariant density n can be eliminated from Hamiltonian by 
solving the second equation and the Hamiltonian becomes a function of p and w. However, 
the constraint is in general unsolvable algebraically. This is the hard part of this issue. 

Now we introduce the major physical assumptions of this paper, under which the existence 
of color skyrmions are considered. 



1. We only consider the fluid configurations with a single characteristic scale, which is 
written as R. 

2. p scales like l/R 3 . 

3. oj scales like l/R. 

4. F follows the 7-law. 



We will simply call R as the size of fluid configuration. Remarkably, these assumptions are 
general enough regardless of the spin and the color of the fluid configuration. Accordingly, 
we introduce the dimensionless quantities in the following: 

X = x/R, P = R 3 P, n = Rto. (17) 
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Moreover, we parameterize the dimensional constant in the definition of 7-law as 

M = aV(4-37). (18) 

So we have two scales in the problem: \x which is fixed, and R which can be varied to minimize 
the energy. It is natural to make the size of the configuration dimensionless by defining 

r = nR. (19) 

Instead of using the invariant density n, it is more convenient to use the ratio 

Z = £. (20) 
n 

With these technical preparations and by putting the 7-law in, we can recast the physical 
quantities in a more compact form: 

H = T + P, 
T = _q_ I d 3 X 1 



a f 

r 5-3 7 J ' 



P 



-r f d 3 X—; (21) 
Z 2 = l + cZ 2 ^" 1 ). (22) 



where we introduce the following definitions 

r 2(4-3 7 ) 



c=-^T (23) 



and 

A = tfP'-r, B=^, C = ^^. (24) 

We will refer T as tension term and P as potential term. We can now define the basic 
problem of the variational procedure: Given 7, whether the total energy in Eq. (21) is able 
to be minimized by the variation of r with the constraint given in Eq. (22)? The solution 
to this question gives rise to a stable configuration in nonabelian fluid within the variational 
method. 

To solve this question, we first consider some general implications of the constraint equa- 
tion (22) and the solution Z as a function in c. Eq. (22) cannot be given algebraically for 
the generic value of 7. In other words, Eq. (22) is transcendental in general; even for the 
simple rational powers like 7 = 7/6, 6/5, 7/5, 8/5, 9/5, 11/6, the roots of this equation cannot 
be written down. Now we rewrite Eq. (22) in another form 

c =f^T ^ 
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we see immediately that Z > 1 for all c > since Z as a ratio of two densities is always 
nonnegative. 1 The next step is to consider the monotonicity of Z in c. We take derivatives 
on both sides of Eq. (25) to get 

dc_ = 2(2- 7 ) 2( 7 - 1) 

dZ Z 2 7"3 ^ Z 2 7-! ' 1 ' 

So as long as 1 < 7 < 2, dc/dZ > and Z is strictly monotonic in c. Henceforth, we will 
assume 

1 < 7 < 2. (27) 
Furthermore, we claim a very useful inequality: 

Z > (7- l)cZ 2 ">- 3 . (28) 
In fact, take derivatives on the both side of Eq. (22) and one can derive the following relation: 

7 2-y-2 

Z _ ( 7 _ 1)cZ 2 7 -3 = __ (29) 

The claimed statement follows the fact that the righthand side of above equation is positive. 
So much for our general consideration of Z. 



3.2 7 = 4/3 

In the following we will discuss the existence of variational solution for different 7. Recall 

the definition c = Cr 2 ^ 7-4 -*. An important observation is that 7 = 4/3 is a critical value 

for which the density ratio Z is independent of the configuration size R. So even though the 

constraint in Eq. (22) can be solved algebraically in this case, the detail of the solution is 

irrelevant to the existence of stable configuration. In fact, the total energy is of the form 

_ const. , . 

E = — — (30) 

where the unspecified constant comes from the dimensionless integral. Physically, this sim- 
plicity follows directly from the fact that the case 7 = 4/3 describes the massless particle 
system. The immediate inference is that no stable fluid configuration with fixed topological 
charge exists at finite R. For the physics of QGP, this is not bad news because we expect to 
see the expansion process of the fireball of nuclear matter. So we refer this type of unstable 
configuration as the unstable color skyrmion. We can imagine that once an unstable config- 
uration is generated, it can be stabilized through a process in which the size R is enlarged 
and the energy is dissipated. Because the process is continuous, the topological charge is not 
changed during this process. However, the final state is a null state in which the configuration 
is completely diluted and the energy is completely dissipated. 



technically, p can be zero at some points. However, such points have no contribution to energy; therefore 
in the integral they can be excluded. 
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3.3 Asymptotic Analysis 



We discuss the asymptotic behaviors of the total energy in the limits r — > and oo for 
1 < 7 < 4/3 and 4/3 < 7 < 2. To do so, we need to solve Eq. (22) approximately in these 
two limits. In fact if c = the physical solution to Eq. (22) is Z = 1. So, for small c, the 
equation can be solved approximately by 



+ (31) 



For large c, we recast Eq. (22) in the form 



= l + Z 2 ^~ 1 \ (32) 



c c 

Let Z ~ c A with a positive power A such that the first term on the right hand side of (32) 
can be dropped. The approximative solution is 

Z = c V2(2-7). (33) 

In the region 4/3 < 7 < 2 for small r, the contribution of Z can be omitted and the 
potential term is the dominant part of the energy in Eq. (21), 



In the region 4/3 < 7 < 2 for large r, 



Z „ r (37-4)/(2- 7 ) p ^ M (35) 

therefore, the tension term dominates. In the region 1 < 7 < 4/3 for large r, the contribution 
of Z can be omitted and again the potential term is dominant. In the region 1 < 7 < 4/3 
for small r from the same equation (35), the tension term dominates. The common feature 
of both regions of 7 is that the energy blows up in the r — > limit and vanishes as r — ► 00. 
So the most immediate conclusion from above asymptotic analysis is that there always exists 
unstable color skyrmion in each region for 7. 

Nevertheless, there is a subtle sub-region for both 1 < 7 < 4/3 and 4/3 < 7 < 2 from 
the behavior of energy in the small r limit. In the sub-region 5/3 < 7 < 2 for small r, the 
tension term behaves like 

T ~ nr^~ 5 0. (36) 

Combining the second formula in Eq. (35) and Eq. (36), we conclude the profile of tension 
term has maximum for some finite r in the region 5/3<7<2. In the sub-region 1 < 7 < 6/5 
for small r, the potential term behaves like 



Mr (6-5 7 )/(2- 7)j (3?) 
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which implies that the profile of potential term has maximum for some finite r in the region 
1 < 7 < 6/5. These results reveal the possibility that there may exist metastable color 
skyrmions in these two sub-regions. By a metastable soliton, we mean that the energy has 
a local minimum. The implication for physics is the following. At the classical level, if the 
configuration is generated with a large enough size R then it is unstable and it can only 
be stabilized by expansion and dissipation until the configuration is completely diluted. For 
small sized configuration however, there is a locally stable point with finite energy. 



3.4 6/5 < 7 < 5/3 

Beyond the asymptotic analysis, we conduct a monotonicity analysis to scrutinize a claim 
that there does not exist metastable color skyrmion in 6/5 < 7 < 5/3. This conclusion is 
certainly true for 7 = 4/3 from our previous discussion. For 4/3 < 7 < 5/3, the conclusion 
is also easy to make. In fact, from the monotonicity of Z in c and the dependency of c on r, 
it is easy to conclude that the total energy is strictly decreasing in r. And the two limiting 
cases r — > 0, 00 have been considered in the asymptotic analysis. 

The case for 6/5 < 7 < 4/3 is more subtle. The monotonicity of total energy is determined 
by the two combinations, 

W 1 = ^^ r , W 2 = f 37 - , (38) 

c 2 7 (4-3 7 ) c 2(4-3 7 )(2- 7 ) 

which are implicit in the expressions of potential term and tension term in total energy. By 
Eq. (22), we know Wi satisfies the relation: 

Wl = C~ 0-r-4)(2-7) + C3^W 2 2(7_1) . (39) 

Taking derivatives on both sides, 

2(W 2 - X^V-) _ (-)- (40) 

v c 4 - 3 i ' ac (4 - 37)c< 4 - 3 ^)( 2 -^) 

where we have expressed the right hand side back in terms of Z. The bracketed factor on 
the left hand side can be also expressed in terms of Z as 

C 4 ~ 3 T c 2(3 7 -4)(2- 7 ) 

The last inequality is because of (28) so we conclude that 
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namely W2 is strictly decreasing in c, hence increasing in r. Due to Eq. (22), W\ satisfies 

3(7-1) 6-57 0^0 
W{ = c 7(4-37) + c 7(4-37) wp . (43) 

Taking derivatives on both sides of Eq. (38) , 

7-1 ^ , dw, (5 7 - e)^- 1 ) + 5i2=il 

2(Wi - ^^^i 7 " 3 )-^ = (")^-^ 3(1-!) C ■ ( 44 ) 

C 7(4-3 7 ) « C _ 3 7 ) C 7(4-3 7 ) 

Again, the sign of dW\/dc is determined by the righthand side because the inequality (28) 
implies the left hand side is positive. So Eq. (44) implies 

namely W\ is also strictly decreasing in c. 2 To conclude this part of the monotonicity 
analysis, we see the total energy is strictly decreasing in r in the region 6/5 < 7 < 5/3. In 
addition, the monotonicity cannot be established from the righthand sides of Eqs. (40,44) 
if 7 > 5/3 or 7 < 6/5. We leave the issue for possible metastable configurations to further 
work. 



3.5 7 = 2 

Now we consider the boundary value of 7 = 2. This is the case which has been worked out 
with a more specific ansatz in our previous paper [8]. The same result can be established on 
more general grounds here. For 7 = 2, Eq. (22) can be solved by 



• I' 1 4- (46) 



The Hamiltonian density in Eq. (14) is then 
The total energy is given by 

H = %J t Px(rtf + ^). (48) 

The existence of color skyrmion that minimizes the energy at certain finite R is straightfor- 
ward. 



2 Certain uniformity conditions for the profile functions Q and P are expected to guarantee that the con- 
clusion does not change after the integral over the space coordinates J d 3 X. 
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3.6 7 = 1 



In this case, Eq. (22) is solved by 




(49) 



and the energy is 




The unstable color skyrmion has a finite energy in the limit R — ► oo and the final state energy 
is 

E — ► jU J d 3 X P. (51) 

To this point, we establish the results of the time- independent color skyrmions listed in 
the introductory section. 



4 7-Dependence of Total Energy 



In previous section, we solve the variational problem for the nonabelian fluid system with 
fixed 7. However, in reality the power 7 is a changing quantity during the process of expansion 
of QGP. So in this section, we consider the dependence of the total energy on the power 7. To 
do so, we will calculate dH/dj. We assume the profiles of £1, P and the parameters //, R are 
independent to 7. We claim that dH/dj < in the "physical parameter region" 1 < 7 < 4/3 
provided the soliton size R is large enough. In fact, 



1 dH 



~ = d>X(( 



„5-3 7 2-2-7 



P7 



(2"7)-sz 



n 2 P 2 -^ 



+ 



p7 



1 In ■ 



5-3 7 ^3-7 r 3( 7 _l)2- 7 +l 



P~ 

dz_ 



p7 



7 2 r 3( 7 -i)z 7 y 

We can impose the following sufficient conditions to guarantee dH/dj < 0: 



^2(7-1) 
r 2(4-3 7 ) 

dz_ 

r 3 Z 



< 



> 



> 



p2( T -l) 



0; 
P. 



(52) 



(53) 

(54) 
(55) 
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By using Eq. (22) and the definition of C in Eq. (24), one can show that (53) is equivalent 

to 



Z 2 < 2+1. (56) 

7 

By the fact that dZ/dc > and the relation in Eq. (23), we conclude that for large enough 
r, the inequality (56), hence (53), is satisfied if 1 < 7 < 4/3. For (54), we take derivatives 
with respect to 7 on both sides of (22), 

(Z - (7 - l)cZ^) ™ = c2F-* In ^. (57) 

Recall the relation in (28). So (54) is satisfied if (55) is satisfied. To show (55), we need to 
check the monotonicity of the combination 

W 3 = c 3/(67-8) Z _ ( 58 ) 

We claim that 

In fact, W3 satisfies the equation 

W 2 = c 3/3 7 -4 + ^37-4^2(7-1) _ (6Q) 

Take derivatives in c: 

2(^3 - ( 7 - l )c 2/37-4w27- 3) ^3 = _ 1_ (3c 7-3 7 /3 7 -4 + 2^-37/37-4^-2). ^ 

dc 4 — 37 

It is easy to show 

W 3 - (7 - l)c 2/37 - 4 W 3 27 " 3 = c 3/67 ~ 8 (Z - (7 - l)cZ 27 " 3 ). (62) 

Therefore, this factor is positive. So from (61), it is obvious that dW^/dc < for 7 < 4/3. 
We see that the conditions (54) and (55) are universally true provided 7 is in the physical 
region. But (53) requires a large configuration size. 

Now we want to estimate the lower bound for r. From (22) and (55), we derive an 
interesting relation: 

Z 2 > 1 + -5-. (63) 

Combined with (56), we have 

1 + ^ < z 2 < 1±1 r 2 > 7O 2 . (64) 

So given an ansatz profile function CI, we know how to estimate the configuration size such 
that the energy is strictly decreasing with 7. 
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The physical meaning of the condition dH/d^ < is the following. It is easy to understand 
from the previous section that with expanding soliton size r, the energy decreases. On the 
other hand, the equation of state in Eq. (2) implies with the decreasing energy density, 7 
becomes smaller as well. So there is a trajectory in the (7, r)-plane along which the total 
energy can be expected to be unchanged during the expansion process. 



5 Time Evolution of Color Skyrmions 

So far we have dealt with time-independent fluid configurations. In this section, we consider 
the time-dependent configurations. To do so, we need to take account of the presence of ujq 
which plays no role in the time-independent case. Accordingly, we need to eliminate j° as 
well as j by the equation of motion, 

5C 



UJ, 



- an*- 2 = (65) 



6j* ^ 

where we have used the 7-law. Eliminating by j^co^ = an 1 ', uj^lo^ = cm 7-2 ^^, we get 

w"^ = aV 7 " 2 . (66) 
The Lagrangian density is therefore expressed in terms of as 

C = ^-an^ = izia-^r^) 3 ^. (67) 

Since in this paper we will concentrate on classical configurations, the pre-factor (7— 1) /^a 1 / 7-1 
does not matter to us at this level. So we will deal with the following Lagrangian 



L 



J d 3 x(u 2 -u}.uj)^. (68) 



As in the case for the time-independent configurations, we make the following scaling 
assumptions: 

where /1 2 are two dimensionless functions depending only on X = x/R. Then the Lagrangian 
in Eq. (68) is transformed to be 

L(R,R) = R^ J d 3 x( fl R 2 -f 2 )^. (70) 

Next we will derive the Euler-Lagrangian equation for the Lagrangian in Eq. (70). First of 
all, 

^ = -L^m k fsx h { h ff-,^, (7i) 
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from which 



-I I d*Xfl{hR 2 -f 2 Y" 

dtdR 7-1^-1 J V ' 



27 ~^R 2 I <?X fl ( fl R 2 -h)^). (72) 



Further 



7 

dL 2 7 -3 n 7^i 



OR 7 

Thus the Euler-Lagrangian equation is given by 

4-37 



Jd?x( h R*-f 2 y^\ (73) 



2 - 

7 



+RRjd*Xf 1 (f 1 tf-f 2 )^ 
+ 2 ^R 2 j d*X h ( h i?- h )^ 

-*LzlJ t p X (f 1 &-f 2 )*^ = 0. (74) 
In the following we will investigate into two particular cases for 7 = 2 and 7 = 4/3. 

5.1 7 = 2 

In this case, stable time-independent color skyrmion exists. So we expect the solution to the 
time-evolution equation to be of the form of oscillations around the stable point. In fact, the 
Euler-Lagrangian equation (74) reduces to 

R 2 

RR + — + 11 = (75) 

where 

fd 3 Xf 2 

" = 57*xa- (76) 

To solve Eq. (75), we change its form to 



3// 



Let q = -R 3 / 2 and p = dq/dt then 
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which can be integrated as 

£ + (79) 

where £ is an integral constant. The solution in Eq. (79) means the color skyrmion for 7 = 2 
forms a one-dimensional Hamiltonian system with the potential of the form q 2 / 3 and the 
motion is always bounded and oscillating! 



5.2 7 = 4/3 

We know that only unstable color skyrmions exist in the time-independent case. Therefore, 
for the time-dependent case, we expect the solution to describe the expansion of the color 
skrymion with the size going from some initial value to infinity. In fact, the Euler-Lagrangian 
equation (74) becomes 

r>2 ±>2 

R 2 (RR- —)- (3 2 {RR- — )+/3 4 = (80) 

where 



" 2 _ o r j% v s2 ' P 4 - 10 f flyrt ' i 81 ) 



fd 3 X fl f 2 fJXf* 
3f<PXf*' P4 12 j(PXf 2 

We will solve Eq. (80) approximately. To do this, we first rewrite this equation by the form 



RR _ 1 = (h ( RR _K_P± (82) 

R 2 4 R2 [ R2 2> R 4 1 j 

and make the following scaling assumptions: 

\R\ > 1, \RR/R 2 \ ~ 1, /? 2 ,4 ~ 1. (83) 

To the zero order R w Eq. (82) can be approximated by 

R<»)R<»> 1 

which is solved by 

R(0) = R^^l + 1)4/3 (85) 

r 

with Rq, t two integral constants. The scaling rule in Eq. (83) works if the ratio 

e = « 1. (86) 

It is very interesting to see that i? 3 / 4 is linear in time. Eq. (85) describes the expansion of 

the color skyrmion for r > and the contraction or collapse for r < 0. Actually, we will 
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only take the former case and consider the latter as unphysical. We will refer the motion in 
Eq. (85) as linear expansion. 

To check the stability of the above-mentioned expansion motion, we do the perturbation 
to the order of e 2 by writing 

R = R(°) + SR (87) 
and introducing the following quantities: 

5R = R e 2 y, x = - + 1, /' = f-. (88) 
r ax 

The perturbation satisfies the following equation 

V»-^ + i»=— (89) 
V 3x + 9x 2 4x 4 /3- l8yj 



Eq. (89) has the particular solution 

y=^l\ (90) 



It is easy to see that the linear independent solutions to the homogenous equation in (89) are 
x 4 / 3 , x 1 / 3 . So the general solution is given by 

y = ^l x m + Ax 4 ^ + Bx 1 / 3 (91) 
8 

where A, B are two integration constants. By examining the initial conditions, we can fix A, 
B. Actually we let 

R\ t=0 = R {0) \t= , R\ t=0 = R^\ t=0 (92) 

then 

y\ x =i, y'\x=i = 0, (93) 

from which we get 

A = ~P2, B = --^ 2 . (94) 
The solution to Eq.(80), to e 2 order, is then given by 

R = R ((l- \fre 2 )x^ + ^/3 2 6 2 x 2 / 3 - \^x^) (95) 

where x = 1 + t/r. To see the stability of the expansion, we consider the long time behavior. 
It is easy to see that the leading order contribution in the limit t — > oo is given by the x 4 / 3 
term, which is the linear expansion as in Eq. (85). 
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6 Conclusion and Discussion 



With the assumption of 7-law for the equation of state for a nonabelian fluid, we show in 
this paper the value of 7 is crucial to the existence and the properties of color skyrmion for 
time-independent configurations. For 7 between 6/5 and 5/3, there is only an unstable color 
skyrmion. For 7 = 2, the skyrmion is stable. The case 7 = 1 is special for the unstable 
skyrmion, since it has finite energy even after it is diluted infinitely. For 1 < 7 < 6/5 and 
5/3 < 7 < 2, we cannot rule out the possibility of the existence of metastable skymions 
besides the unstable ones. For two particular values 7 = 4/3,2, we furthermore consider the 
time-dependent configurations. And for the latter we find the oscillating evolution and for 
the former we find linear expansion. 

We now turn to the question of how this relates to QGP. As mentioned in the introduction, 
if we try to model the lattice estimate of the equation of state as a 7-law, the value of 7 varies 
between 4/3 and 1, as the energy density decreases. Therefore, our analysis of the expanding 
soliton for 7 = 4/3 takes on a special significance. In the creation of the QGP by nuclear 
collision, we are starting at high energy densities, or at values of 7 close to 4/3. The collision 
process also favors the creation of small solitons during the phase transition process since 
establishing coherent fields or color densities over large scales is generally more difficult, 
less probable. Based on these two observations and in the light of our general analysis, an 
expanding soliton would seem to be the generic case we can expect for the QGP. There 
are, of course, some caveats to our analysis. We have neglected thermal gradients as well 
as dissipative effects such as viscosity. (The latter quantity, although small, is not zero.) 
It would be interesting to incorporate some of these effects. The soliton by virtue of its 
nontrivial topology is a generic feature of nonabelian fluid dynamics. We may expect many 
qualitative aspects to hold true even with thermal gradients, viscosity, etc., although there 
can be significant changes in details. 

We note that the possibility of color skyrmions arises also in the recent effective action 
description of Wilson lines [10, 11]. It is clear that skyrmions in the QGP merit further 
analysis. Other hydrodynamics approach towards QGP can be found for example in [12]. 
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